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Long QT syndrome (LQTS) is a congenital abnormality of cardiac repolarization that 
manifests as a prolonged QT interval on 12-lead electrocardiograms (ECGs). The syndrome 
may lead to syncope and sudden death from ventricular tachyarrhythmias known as 
torsades de pointes. An increased persistent Na + current is known to cause a Ca 2+ 
overload in case of ischemia for example. Such increased Na + persistent current is also 
usually associated to the LQT3 syndrome. The purpose of this study was to investigate the 
pathological consequences of a novel mutation in a family affected by LQTS. The impact 
of biophysical defects on cellular homeostasis are also investigated. Genomic DNA was 
extracted from blood samples, and a combination of PCR and DNA sequencing of several 
LQTS-linked genes was used to identify mutations. The mutation was reproduced in vitro 
and was characterized using the patch clamp technique and in silico quantitative analysis. 
A novel mutation (Q1476R) was identified on the SCN5A gene encoding the cardiac Na + 
channel. Cells expressing the Q1476R mutation exhibited biophysical alterations, including 
a shift of SS inactivation and a significant increase in the persistent Na + current. The 
in silico analysis confirmed the arrhythmogenic character of the Q1476R mutation. It 
further revealed that the increase in persistent Na + current causes a frequency-dependent 
Na + overload in cardiomyocytes co-expressing WT and mutant Na v 1.5 channels that, in 
turn, exerts a moderating effect on the lengthening of the action potential (AP) duration 
caused by the mutation. The Q1476R mutation in SCN5A results in a three-fold increase 
in the window current and a persistent inward Na + current. These biophysical defects 
may expose the carrier of the mutation to arrhythmias that occur preferentially in the 
patient at rest or during tachycardia. However, the Na + overload counterbalances the 
gain-of-function of the mutation and is beneficial in that it prevents severe arrhythmias 
at intermediate heart rates. 
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INTRODUCTION 

In cardiomyocytes, action potentials (AP) ensure the propagation 
of excitation and the triggering of contraction. AP characteristics 
are directly related to the activity of voltage-gated sodium (Na + ), 
calcium (Ca 2+ ), and potassium (K + ) channels (Amin et al., 
2010). After a slight depolarization, voltage-gated Na + channels 
(Na v 1.5), encoded by the SCN5A gene, activate, which leads to an 
inward Na + flow followed by a fast depolarization (Gellens et al., 
1992). The AP is composed of a maintained depolarization, the 
"plateau" phase, followed by a final repolarization. The depolar- 
ization at the plateau is supported by a slowly decaying inward 
calcium current, background inward currents, and a persistent 
sodium current. The plateau is also supported by a delay in the 
onset of repolarizing potassium currents. The late activation of 
these potassium currents ensures final repolarization and causes 



the membrane potential to return to the diastolic level after about 
400 ms. The maintenance of action potential duration (APD) is 
critically related to the balance between depolarizing currents 
and repolarizing currents, so that an increase in depolarizing 
current(s) or a decrease in repolarizing current(s) will increase 
APD and, consequently, the QT interval of the electrocardiogram 
(ECG). 

Long QT syndrome (LQTS) includes a congenital abnor- 
mality of cardiac repolarization that causes syncope and sud- 
den death (Huang et al., 2011). This hereditary cardiac disease 
manifests as an increase in the QTc value over 450 ms on 12- 
lead ECGs, which points to an increase in the duration of 
ventricular APs. Under optimal conditions, QTc prolongation 
can induce ventricular tachyarrhythmia called torsade de pointes 
(TdP), which can degenerate into lethal ventricular fibrillation. 
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To date, mutations in thirteen genes have been identified as 
causing 13 types of LQTS (LQT1-13) (Schwartz et al, 2012) 
and are divided into three major groups: (1) mutations in 
genes encoding voltage-gated channels [KCNQ1 ( LQT 1 ), KCNH2 
(LQT2), SCN5A (LQT3), KCNJ2 (LQT7), KCNJ5 (LQT13)], 
(2) mutations in genes encoding subunits regulating voltage- 
gated channel [KCNE1 (LQT5), KCNE2 (LQT6), SCN4B (LQT 
10)], and (3) mutations in genes encoding channel interact- 
ing proteins [ANKB (LQT4), KCNE1 (LQT5), KCNE2 (LQT6), 
CAV3 (LQT9)]. 

SCN5A gene mutations can lead to several inherited cardiac 
diseases such as LQT3, Brugada syndrome (BrS) (Moreau et al., 
2012), progressive cardiac conduction defect (CCD), sinus node 
disease, and dilated cardiomyopathy (DCM) (Amin et al., 2010; 
Gosselin-Badaroudine et al, 2012). In LQT3, a gain-of-function 
mediated by the appearance of a persistent current usually char- 
acterizes Na + channel mutations. The persistent Na + current has 
previously been shown to induce Ca 2+ overload (Tang et al., 
2012). Indeed, in hypoxia condition, the persistent Na + cur- 
rent is increased and have been identified as a predominant 
cause of the deleterous Ca 2+ overload (Tang et al., 2012). Such 
phenomenon could be important in the development of LQT3 
syndrome and would be driven by the activation of Na + /H + 
and Na + /Ca 2+ exchangers. These exchangers are predominantly 
implicated in the Ca 2+ , Na + , and pH regulation. However, 
most mutations in K + channels or their regulatory subunits 
cause a loss-of-function characterized by a dominant negative 
effect. These biophysical defects lead to AP prolongation, creat- 
ing a substrate for the development of ventricular arrhythmias 
(Schwartz et al, 2012). 

The purpose of our study was to investigate the impact of a 
mutation causing LQTS in a patient and her family members 
using an in vitro heterologous expression system, the voltage- 
clamp technique, and in silico analysis. The biophysical investi- 
gation of the mutation in tsA201 cells revealed a positive shift in 
steady-state (SS) inactivation that results in an increased window 
current associated with the presence of a persistent Na + current, 
which would explain the link between the mutation and LQTS in 
the family members. Furthermore, in silico experiments revealed 
that the persistent Na + current induces an increase in the intracel- 
lular Na + concentration. This causes an increase in repolarizing 
currents that shortens APD, counterbalancing the direct effect 
of the gain-of-function caused by the mutation and potentially 
benefitting the patient. 

METHODS 
CLINICAL EVALUATION 

The proband and her family members underwent a detailed clin- 
ical and cardiovascular examination, including a 12-lead ECG, 
a determination of hemodynamic parameters, and laboratory 
blood tests. Echocardiography was used to exclude any structural 
heart disease. The QT interval was measured in lead II or V5 of 
the ECG and was corrected for the heart rate (QTc) using Bazett's 
formula. Exercise testing was performed to detect QT adaptation 
to heart rate changes, which is known to predict the genotype 
and direct therapy in LQTS (Schwartz et al, 2001; Wong et al., 
2010). 



The investigation conformed to the principles set out in the 
Declaration of Helsinki and was approved by the local institu- 
tional ethics committees on human subject research. The partici- 
pants provided written informed consent prior to the study. 

GENOTYPING 

The samples were sent to Genaissance Familion® (New Haven, 
CT, USA) for genetic testing of inherited cardiac syndromes. 
Genomic DNA was extracted and was amplified by PCR. 
Templates of the target exons, splice junctions, and flanking 
regions of KCNQ1 (LQT1), KCNH2 (LQT2), SCN5A (LQT3), 
KCNE1 (LQT5), KCNE2 (LQT6), CACNA1C (LQT8, exons 8 and 
9), CAV3 (LQT 9), SCN4B (LQT 10), AKAP9 (LQT 11, exon 18), 
and SNT12 (LQT 12) were generated and genotyped. 

MUTAGENESIS 

Mutant Na v 1.5/Q1476R was generated using QuickChange 
TM site-directed mutagenesis kits according to the manu- 
facturer's instructions (Stratagene, La Jolla, CA, USA). The 
hNa v 1.5/mutants were constructed from WT channels (NCBI 
accession number: M77235) using the following mutagenic sense 
and antisense primers (mutated sites are underlined): 

1. 5'-CAACTTCAACCAACGGAAGAAAAAGTTAGGG-3' 

2. 5'-CCCTAACTTTTTCTTCCGTTGGTTGAAGTTG-3' 



Mutant and the wild-type (WT) Na v 1.5 in a pcDNA3 construct 
were purified using Qiagen columns (Qiagen Inc., Chatsworth, 
CA, USA). 

TRANSFECTI0N OF THE tsA201 CELL LINE 

TsA201 is a mammalian cell line derived from human embry- 
onic kidney HEK 293 cells by stable transfection with SV40 
large T antigen (Margolskee et al., 1993). The tsA201 cells were 
grown in high glucose DMEM supplemented with fetal bovine 
serum (10%), L-glutamine (2mmol/L), penicillin (100 U/ml), 
and streptomycin (lOmg/ml) (Gibco BRL Life Technologies, 
Burlington, ON, Canada) and were incubated in a 5% CO2 
humidified atmosphere. The cells were transfected using the 
calcium phosphate method (Margolskee et al, 1993), with the 
following modification to facilitate the identification of indi- 
vidual transfected cells: the cells were co-transfected with an 
expression plasmid for a lymphocyte surface antigen (CD8-a) 
(Jurman et al., 1994). The human Na + channel Pi subunit and 
CD 8 were constructed in the piRES vector (piRES/CD8/Pi ). With 
this strategy, transfected cells that bind beads also express the 
Pi -subunit. cDNA (5u,g) coding for the WT or mutant Na + 
channel a-subunit and 5 |ig of piRES/CD8/Pi were used. For 
the patch clamp experiments, 2 to 3-day-post-transfection cells 
were incubated for 5min in a medium containing anti-CD8-a 
coated beads (Jurman et al, 1994) (Dynabeads M-450 CD8-a). 
Unattached beads were removed by washing. The beads were pre- 
pared according to the manufacturer's instructions (Dynal A.S., 
Oslo, Norway). Cells expressing surface CD8-a fixed the beads 
and were visually distinguishable from non-transfected cells by 
light microscopy. 
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PATCH CLAMP TECHNIQUE 

Macroscopic Na + currents from tsA201-transfected cells were 
recorded using the whole-cell configuration of the patch clamp 
technique (Hamill et al., 1981). Patch electrodes were made from 
8161 Corning borosilicate glass and were coated with HIPEC 
(Dow-Corning, Midland, MI, USA) to minimize their capaci- 
tance. Patch-clamp recordings were made using low resistance 
electrodes (<1M£2), and a routine series resistance compensa- 
tion by an Axopatch 200 amplifier (Axon Instruments, Foster 
City, CA, USA) was performed to values >80% to minimize 
voltage-clamp errors. Voltage-clamp command pulses were gen- 
erated by microcomputer using pCLAMP software vlO.O (Axon 
Instruments). Na + currents were filtered at 5 kHz, digitized at 
10 kHz, and stored on a microcomputer equipped with an AD 
converter (Digidata 1300, Axon Instruments). 

SOLUTIONS AND REAGENTS 

For whole cell recordings, the patch pipette contained (in 
mmol/L): 35 NaCl, 105 CsF, 10 EGTA, and 10 Cs-HEPES. The pH 
was adjusted to 7.4 using 1 N CsOH. The bath solution contained 
(in mmol/L) 150 NaCl, 2 KC1, 1.5 CaCl 2 , 1 MgCl 2 , 10 glucose, and 
10 Na-HEPES. The pH was adjusted to 7.4 with 1 N NaOH. The 
recordings were made 5 min after establishing the whole cell con- 
figuration to allow the current to stabilize and to permit adequate 
diffusion of the contents of the patch electrode. All the recordings 
were made in the following order: I/V curve, SS inactivation, and 
recovery from inactivation. Experiments were carried out at room 
temperature (22-23°C). 

In silico analysis and AP modeling 

The Priebe and Beuckelmann (1998) (PB), the Iyer et al. (2004) 
(IMW) model and the Ten Tusscher et al. (2004) (TNNP) mod- 
els were considered as possible alternatives to represent human 
ventricular myocyte (Priebe and Beuckelmann, 1998; Iyer et al., 
2004; Ten Tusscher et al, 2004). A comparison of these models 
was done by Ten Tusscher et al. in 2006 and it showed that the 
TNNP model was the most pertinent to reproduce the effects of 
IKr and IKs block on the human ventricular AP duration (Ten 
Tusscher et al, 2006). We considered it as a valid model to address 
changes related to LQT conditions. Besides, rather than using 
Markov models, this model is based on Hodgkin-Huxley equa- 
tions, which are easy to adjust to known experimental data (Fink 
and Noble, 2009). 

Thus, the in silico analysis was performed with the Ten 
Tusscher et al. (2004) model (Ten Tusscher et al., 2004), using the 
tentusscher_noble_noble_panfilov_2004_a.cellml file, ran in the 
COR environment (Peter Garny, http://cor.physiol.ox.ac.uk/). 

Changes from the original set of equations are described in the 
Appendix. 

When we attempted to modify the activation and inactiva- 
tion parameters of the Na + current (INa) in the original TNNP 
model, instabilities were found, which could not be resolved. 
Thus, the mathematical equations used for activation, inactiva- 
tion, and recovery from inactivation of INa were taken from a 
stable guinea-pig ventricular cell model adjusted for 37°C (Pasek 
et al., 2008). However, the processes were accelerated two-fold 
(the time constants were reduced two-fold). This was done to 



allow the model to reproduce the experimentally observed Na + 
currents when we ran, in the model set for room temperature, the 
same voltage-clamp protocols as used in our in vitro experiments. 

To adapt to the formalism of the Na + current (INa) calcu- 
lation that uses m 3 , the measured values of the SS apparent 
activation variable measured in tsA201 cells were fitted with the 
following equation: m m f = [1/(1 + exp((V m — Vh)/fc))] 3 , which 
yielded values of Vh = —59.1 mV and k = —7.9 mV for WT and 
V h = -56.8 mV and k = -8.7 mV for Q1476R. This explains 
why these values of Vh and k were different from those calcu- 
lated by fitting a simple Boltzmann function to data measured in 
tsA201 cells (Table 1). 

Lastly, 5.5 mV was added to Vh and 0.5 mV to k to account 
for the difference in temperature given that the experiment was 
conducted at room temperature whereas the model was designed 
toworkat33°C, such that V h = -59.1 + 5.5andfc= -7.9 + 0.5 
for WT and V h = -56.8 + 5.5 and k = -8.7 + 0.5 for Q1476R. 

The SS inactivation functions were h^r = 1/(1 + exp((V m — 
V h ) A)), where V h = -92.0 + 5.5 and k = 4.4 for WT, and V h = 
-86.1 + 5.5 and k = 4.3 for Q1476R. In addition, 5.5 mV was 
added to Vh to account for the difference in temperature. 

This adjustment of the kinetic parameters was derived from 
the detailed analysis of the effects of temperature on the activation 
and inactivation kinetics of Na v 1.5 channels by Nagatomo et al. 
(1998). It reflects a change in temperature from 23 to 33°C. A lin- 
ear extrapolation to 37°C could have been done. Nevertheless, no 
differences should be observed as these corrections would apply 
to both WT and mutated fast INa model equations. 

The modeling of the heterozygous condition is based on the 
assumption that, in the carrier, both the WT and the mutated 
genes would lead to the presence of equal amounts of each cor- 
responding channel protein. Thus, half of the fast Na + current 
is characterized by the parameters of the WT channels while the 
second half is characterized by the parameters of the Q1476R 
channels. 

The persistent part of the fast Na + current was initially absent 
in the TNNP model. A persistent current was added that would 
not inactivate but that would activate and deactivate following 



Table 1 | Biophysical properties of Na v 1.5/WT and Q14767R. 







WT(n = 13) 


Q1476R (n= 14) 


Peak current (pA/pF) 




-739.1 ±73.1 


-690.8 ±69.9 


Persistent current (% 


I Na+) 


0.36 ± 0.08 


0.74 ± 0.11* 


Activation 


1/1/2 (mV) 


-47.9 ± 1.2 


-46.3 ± 1.1 




*: (mV) 


-5.6 ± 1.2 


-6.5 ±0.4 


Inactivation 


1/1/2 (mV) 


-92.0 ± 1.3 


-85.5 ±0.4* 




*: (mV) 


4.4 ±0.1 


4.3 ± 0.1 


Recovery from fast 


A 


99.9 ±0.05 


98.7 ± 1.1 


inactivation 


Tau 


2.7 ±0.1 


2.2 ± 0.3 



Vj/2, Mid-point for activation or inactivation. 
k, Slope factor for activation or inactivation. 
'Difference between Na v 1.5/WT and Na v 1.5/Q1476R (p < 0.05). 
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B QTc 

Baseline 
490 msec 



I 



II 




unaffected 
daughter 
425 msec 



FIGURE 1 | Family pedigree and 12-lead ECG of the proband. (A) Family 
pedigree. Subjects II-2 and III-2 and III-3 with the LQT3 phenotype (black 
circle and square), the arrow indicates the proband. (B) Twelve-lead ECG. 
Baseline ECG in the proband, demonstrating a long isoelectric ST segment, 
followed by a deferred T was with normal morphology. The second ECG 
represents a repeat ECG after 2 doses of 200 mg of mexilitine. The ECGs of 
her affected and unaffected daughter are also presented. Voltage scale: 
0.5mV/large square, time scale: 0.2s/large square. 



the same equations as the fast Na + current (see equations in the 
Appendix). 

Initially, the model was run with the basic set of values of vari- 
ables and constants as used in the TNNP 2004a model, except for 
the formulation of INa, as described above. The APD90 was mea- 
sured as the time between upstroke and 90% repolarization from 
the highest level of the plateau. No attempt was made to readjust 
the APD90 at SS at 1 Hz, which was increased from 322 ms in the 
original TNNP model to 454 ms for our WT settings. 

Voltage-clamp testing was performed using MATLAB to show 
that the basic properties of the fast INa current, the persistent Na + 
current and the window current, adjusted for room temperature, 
matched those of the WT and mutated currents as determined 
in the expression experiments when using the same protocols 
and measurement procedures in the presence of a realistic series 
resistance (not shown). 

Based on our previous modeling experience with a persistent 
INa (Christe et al, 2008), it was necessary to allow the model to 
stabilize with a small amount of persistent INa before increasing 
the persistent INa until reaching at most one-third of the frac- 
tion of the fast peak INa observed in expression experiments. This 
was the maximal amount of persistent INa tolerated by the het- 
erozygous mutation model in order to ensure that a stable AP was 
generated at 1 and 1.5 Hz. 

DATA ANALYSIS AND STATISTICS 

The electrophysiological data were analyzed using macros 
in Clampfit (pCLAMP vlO.O, Molecular Devices) and cus- 
tom programs written using MATLAB (The MathWorks Inc.). 
Data are expressed as means ± sem (standard error of the 
mean). Statistical comparison were performed using an unpaired 
Student's f-test or One-Way Anova with Bonferroni post-hoc 
test in Sigma plot (Jandel Scientific software). Differences were 
deemed significant at a P-value < 0.05. 

RESULTS 

CLINICAL CHARACTERISTICS OF THE INDEX PATIENT AND THE FAMILY 
MEMBERS 

A 34-year-old female was assessed, in the context of a family 
history of sudden death, for a remote history of syncope after 
awakening to answer the telephone. Her sister had experienced a 
series of similar episodes at age 14 and had subsequently collapsed 
at rest and could not be resuscitated. The autopsy did not reveal 
a cause of death. A maternal cousin and grandaunt also experi- 
enced sudden death (Figure 1A). The resting QTc of the proband 
ranged from 450 to 485 ms and shortened to 422 ms at peak exer- 
cise. The echocardiogram was normal. Nadolol 40 mg once daily 
and mexiletine 200 mg twice daily were initiated. While mexile- 
tine shortened the QTc (Figure IB), it was discontinued before 
reassessment because of side effects. We identified an A-to-G base 
change at nucleotide position 5300 in exon 25 on the SCN5A 
gene by PCR and DNA sequencing. The base change led to a 
glutamine (Q)-to-arginine (R) substitution at residue 1476. The 
Q1476 residue is in a region of the SCN5A gene/Na v 1.5 pro- 
tein in Na + channels (Figure 2A) that has been extraordinarily 
highly conserved throughout evolution, from humans to zebra 
fish (not shown). Another highly conserved glutamine (Q1475, 



Figure 2A) was located upstream from the Q1476 residue. The 
mutation occurred in an important region of the Na + channel 
known as the inactivation gate (Figure 2B). An H558R polymor- 
phism was associated with the novel Q1476R mutation. The four 
children of the index patient were examined and presented QTc 
values of 414, 444, 480, and 502 ms. Three of them underwent an 
exercise test and their QTc values went from 444 to 442 ms, 480 to 
500 ms, and 502 to 522 ms, respectively. The two who presented 
with borderline QTc prolongation were subsequently confirmed 
to be carrying the Q1476R variant. The index patient remained 
free of symptoms during the 5 -year follow-up period. 

BIOPHYSICAL CHARACTERIZATION 

We investigated whether the novel mutation caused an alter- 
ation of the biophysical properties of the Na v 1.5 Na + chan- 
nel that would explain the clinical phenotype of the index 
patient. Macroscopic Na + currents were recorded from tsA201 
cells expressing either WT (Na v 1.5/WT) or mutant channels 
(Na v 1.5/Q1476R) co-transfected with the regulatory Pi-subunit 
(see Methods for more details) (Figure 3A). The mutant Na + 
channel current exhibited no significant shift in the current- 
voltage (I/V) relationship, and the current density was similar 
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DIIIS6 



DIVS1 
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SCN3A_Human 
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SCN5A Human 

SCN8A_Human 

SCN9A_Human 

SCNIOAHuman 
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LFIGVIIDNFNQQKKKFGGQDIFMTEEQKKYYNAMKKLGSKKPQKPIPRPGNKFQGMVFDFVTRQVFDI 
LFTGVI IDNFNQQKKKFGGQDIFMTEEQKKYYNAMKKLGSKKPQKPIPRPANKFOGMVFDFVTKQVFDI 
LFIGVIIDNFNQQKKKFGGQDIFMTEEQKKYYNAMKKLGSKKPQKPIPRPANKFQ.GMVFDFVTRQVFDI 
LFIGVIIDNFNQQKKKLGGKDIFMIEEQKKYYHAMKKLGSKKPQKPIPRPQNKIQGMVYDLVTKQAFDI 
LFIGVI IDNFNQQKKKLGGQDI FMTEEQKKYYNAMKKLGSKKPQKPIPRPLNKYQGFIFDIVTKQAFDV 
LFIGVI IDNFNOQKKKFGGQDI FMTEEQKKYYNAMKKLGSKKPQKPI PRPLNKIOGIVFDFVTQOAFDI 
LFIGVIIDNFHQQKKKLGGQDIFMTEEQKKYYHAMKKLGSKKPQKPIPRPGHKIQGCIFDLVTHQAFDI 
LFVGVI IDNFHQQKKKLGGQDI FMTEEO.KKYYNAMKKLGSKKPQKPIPRPLNKFQGFVFDIVTRQAFDI 
LFIGVIIDHFHQQQKKLGGQDIFMTEEQKKYYHAMKKLGSKKPQKPIPRPLHKCQGLVFDIVTSQIFDI 




] Voltage Sensitive Domain 



I I Pore Domain 

Inactivation gate containing IFM motif (DIII-DIV linker) 

+ Q1475 (FHM) highly conserved amino acid. 
+ Q1476 (LQT3) highly conserved amino acid. 



FIGURE 2 | Q1476 is within a highly conserved region. (A) Amino 
acid alignments of the region between transmembrane domains IIIS6 
and IVS1 (lll-IV linker) showing in red that glutamine 1475 and 
glutamine 1476 are highly conserved in human voltage-gated sodium 
channels. The IFM inactivation particle is displayed in red. (B) Schematic 
representation of the Na v 1.5 channel structure displaying voltage 



sensitive domains, the segments forming the pore domain, and the 
inactivation gate. Q1475 and Q1476 are also indicated. The missense 
mutation of one these two amino acids leads to several pathologies, 
depending on the channel impacted. The mutation of Q1475 in Na v 1.1 
causes familial hemiplegic migraine (FHM), whereas the mutation of 
Q1476 in Na v 1.5 causes LQT3 syndrome. 



to that of the WT channel (WT: -739.1 ± 73.1 pA/pF, Q1476R: 
-690.8 ± 69.9 pA/pF) (Figure 3B, Table 1). This result is based 
on the measurement of the peak current which is defined as 
the maximal absolute current obtained with the stimulation pro- 
tocol (see Figure 3 legend). This peak current corresponds to 
the value where the global conductance reaches its maximum. 
A small persistent Na + current was recorded. The persistent 
Na + current was investigated using a 400 ms two-pulse pro- 
tocol, and it was measured as mean of the current between 
100 and 350 ms. The magnitude of this persistent current was 
significantly (p < 0.05) different from that recorded from WT 
channels at -30 mV (0.36 ± 0.08%, n = 13 for WT channels, 
0.74±0.11%, n = 14 for mutant channels) (Figure 3C, Table 1). 
SS activation (Gv) and inactivation curves were also studied 
(Figure 3D, Table 1). The Gv curve of the mutant Na + channel 
showed no significant shift (p = 0.340 for V1/2 andp = 0.081 for 
the slope factor). However, there was a 6.5 mV (p < 0.05) shift 
in the SS inactivation curve toward more positive voltages for 
the Q1764R mutant channel (Vy 2 Na v 1.5/Q1476R = -85.5 ± 
0.4 mV, n = 14 vs. V 1/2 Na v 1.5/WT = -92.0 ± 1.3 mV, n = 13). 
The slope factor was not significantly affected (p = 0.544). The 
mutation did not induce any significant changes in the time 
course of the Na + current decay or the recovery from inactivation 
(Figure 4). 

The overlap of SS activation and inactivation of Na + channels 
defines a range of voltages (i.e., window) where Na + channels 
open, resulting in an inward Na + current that could potentially 
depolarize the resting membrane potential and increase myocyte 



excitability. The Q1476R mutation induced an increase in the 
overlap of Na + channel activation and inactivation by shifting 
the inactivation to more positive potentials. This shift tended to 
expand the window and thus, the fraction of permanently acti- 
vated channels. Figure 5A shows the predicted window currents 
of the WT and Q1476R mutant Na v 1.5 channels. The Q1476R 
mutation produced a ~7-fold increase in the Na v 1.5 window 
current, suggesting that more channels had reopened within the 
window region. Figure 5B shows the increase in the window cur- 
rent together with the increase in the perisistent component, 
suggesting that Q1476R had a higher probability to open at 
voltages ranging from — 140 to +20 mV. 

The persistent Na + current can be attributed to the Na v 1.5 
channel as it was blocked (0.06 ± 0.06%, n = 5) by 10 |xM of 
tetrodotoxin, a Na + channel blocker (Figure 6). Mexiletine was 
seen to be efficient to treat the patient. The effect of this class 
I B antiarrhythmic was thus tested on the persistent Na + cur- 
rent. 50 ixM of mexiletine importantly reduced the persistent Na + 
current (0.09 ± 0.03%, n = 5) (Figure 6). 

IN SILICO CHARACTERIZATION OF THE CHANGES INDUCED BY THE 
HETEROZYGOUS MUTATION Q1476R 

The in silico characterization incorporates the heterozygous 
aspect of the mutation. This model is used to observe the effects 
of the mutation on several parameters including AP, Na + window 
current INa (the large inward current surge simultaneous with 
AP upstroke is largely off scale), persistent Na + current (INap), 
Na + and Ca 2+ intracellular concentrations (Na; and Ca;). As a 
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FIGURE 3 | Biophysical characterization of Na v 1.5/WT and the 
Q1476R mutant channels. (A) Representative whole-cell current 
traces of the WT (left) and Q1476R (right) channels. Currents were 
elicited using a voltage-clamp protocol where depolarizing pulses 
were applied for 50 ms from —100 to +80 mV in 5mV increments 
(see protocol in inset). The zero current corresponds to the basal 
line shown at the beginning of each trace. (B) After normalization to 
cell capacitance, current-voltage (l-V) relationships were constructed 
for the Na v 1.5/WT (n=13, black) and Q1476R (n = 14, red) channels. 
(C) Representative current traces obtained after a — 30 mV 
depolarizing pulse showing INap. The Na v 1.5/WT channel trace is 
depicted in black while the Na v 1.5/Q1476R channel trace is depicted 



in red. (D) Voltage dependence of SS activation (circles) or 
inactivation (squares) for the WT (black) and Q1476R (red) Na v 1.5 
Na v 1.5 channels. A standard Boltzmann function was used to fit 
activation curves and to obtain the V-\/2 and given in Table 1 
[G(V)/G max = 1/(1 + exp(— (V — Vi/2)/k))]. Inactivation currents were 
obtained by applying conditioning pre-pulses to membrane potentials 
ranging from a holding potential of —140 to — 10 mV for 500 ms in 
5mV increments and were then measured using a 4-ms pulse to 
— 30 mV at each step (see protocol in inset). The recorded 
inactivation values were fitted to a standard Boltzmann equation 
WVO/'max = 1/(1 + exp((V- V y2 )/k)) + C] to obtain the values given 
in Table 1. 



rule, model simulations were conducted for at least a 2000 s cell 
lifetime to ensure that SS in all variables was reached. 

With 1-Hz stimulations, at SS, the WT model shows an APD90 
of 454 ms and a Nai of 11.94mmol/L (central column, trace a, 
in Figure 7 A, and Table 2). The model was then changed to the 
heterozygous condition (denoted as WR), with half of the INa 
current being WT and half being Q1476R. The first AP obtained 
in the WR condition had an APD90 of 587 ms (trace b) without 
appreciable change in Nai. Both Na + window current and INap 
amplitudes were increased when the mutation was introduced 



(2nd and 3rd panels from the top, respectively, central column 
of Figure 7A). This was consistent with estimates computed from 
the voltage-clamp data (Figure 5). The difference is that the mag- 
nitude of the mutation-induced increase in the window current 
and in the INap as evaluated in the heterozygous model was 
two-fold less than that reported in Table 1, since, in the WR 
model, only half of the current is of the mutated type. When 
running this heterozygous model to SS (trace c), Na; expect- 
edly increased, reaching 12.92 mmol/L, while APD90 shortened 
to 513 ms. The diastolic Ca; and Ca; transients were also increased 
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FIGURE 4 | Time constants of inactivation and of recovery from 
inactivation. (A) Curves representing the time constant of inactivation as a 
function of the voltage for the WT and Q1476R channels. A two-exponential 
function was used to fit the inactivation curves observed on whole cell 
current traces: / = Aast(exp[-(f - fc)/tfastl) + Asiow(exp[-(f - fc)/t S !ow] 
+ 0, where /\ fast and A S | 0W are fractions of recovery of the fast and slow 
components, f is time, and k is the delay factor for activation or inactivation. 
tfast and x s | ow were obtained from the equation from —50 to +20 mV for 
the Na v 1.5/WT and Q1476R channels. (B) Recovery from fast inactivation of 
the WT (black) and Q1476R (red) channels. Curves were obtained using the 
two-pulse protocol shown in the inset. The pulses (+30 mV) induced 
maximal activation, and the time interval between pulses represents the 
time needed to obtain a re-activation. The recorded recovery from 
inactivation values were fitted to a single exponential equation 
[/4*(1 — exp(— x/x) + C)l to obtain the values given in Table 1. 



(Figure 7 A, lower panel, central column and Table 2). The instan- 
taneous APD90 lengthening by 29% due to changing from WT to 
WR condition was reduced to 13% vs. WT when the WR model 
reached SS. Most of the lengthening was thus counterbalanced 
by changes that occurred during regular stimulations for longer 
periods (typically 2000 s). 

The same simulations were conducted at 1.5 Hz until SS 
(Figure 7 A, 3rd column, and Table 2). APD90 and Na;, respec- 
tively stabilized at 392 ms and 13.59 mmol/L for WT (trace a). 
After the switch to the WR model, the first AP (trace b) had 
an APD90 of 471 ms (Naj remained unchanged). After reach- 
ing SS (trace c), APD90 shortened to 436 ms and Na; increased 
to 14.57 mmol/L. At SS, the diastolic Ca; and Ca; transients 
had increased in the WR model. The increases were larger than 



Na v 1.5/Q1476R 
Na„1 .5/WT 




V(mV) 

FIGURE 5 | Na v 1.5/WT and Q1476R window currents. (A) The window 
current is defined as the overlap between activation and inactivation. The 
window current was obtained using the following equation: (1/(1 + 

exp((Vi ^activation - W/CactivatioJ) X ((1 - 0/(1 +exp((V- V, /2 inactivation )/ 

inactivation)) + O. (B) The INap component was added to the window 
current of the WT or the Q1476R channel. 



those seen at 1 Hz (Figure 7A, lower panel, 3rd column and 
Table 2). 

To test the effect of bradycardia, the stimulation rate was 
lowered from 1 Hz to 0.5 Hz. At SS, the WT model had 
a longer APD 90 ( 5 58 ms) and a lower Na; (9.84 mmol/L) 
than at 1 Hz. In the WR model at 0.5 Hz, APD 90 and Na; 
reached SS values of 705 ms and 10.95 mmol/L (Table 3), while 
the diastolic Ca; and Cai transients decreased (Figure 7A, 
1st column, lower panel and Table 3). To determine the 
effect of changes in Nai on APD90, changes in Na; were 
applied instantaneously to the model in several settings 
(Figure 7B). 

In the WT model, from SS at 1 Hz, setting Na; to its value at SS 
at 0.5 Hz produced an APD90 of 542 ms (Figure 7B, left panel). 
This was comparable to the SS APD90 at 0.5 Hz when Nai was 
not imposed (558 ms). Likewise, when SS Nai wa s imposed in 
the WR model at 0.5 and 1 Hz, the APD 90 s of the next AP were 
respectively 683 and 515 ms, vs. 705 and 513 ms at SS when Na; 
was not imposed (left panel of Figure 7B and Table 3). At 1 Hz, 
in WR model under SS, changing Naj to WT SS value produced 
APD90 of 582 ms. This value is close to that of the first WR AP 
(APD90 of 587 ms). Then, changing to WT produced APD90 of 
45 1 ms which is similar to SS WT APD90 of 454 ms (middle panel 
of Figure 7B and Table 1). At 1.5 Hz (Figure 7B, right panel and 
Table 2), APD90 produced by changing Na; values were very close 
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FIGURE 6 | Effect of Na+ channel blockers on the INap. 

(A) Representative current traces obtained after a —30 mV depolarizing 
pulse showing INap before and after drug treatment. The Na v 1.5/Q1476R 
traces are depicted in black while the traces after treatment are shown in 
red for mexiletine (50 \lM, left) and blue for TTX (10 \lM, right). (B) 
Histogram of INap. The INap accounted for 0.36 ± 0.08% (n= 13, red 
column) of the peak current for WT channel and 0.74 ± 0.11 % (n = 14, 
orange column) for Q1476R channels ( # p = 0.015 when compared to WT 
channels). The treatment with 50 |iM of mexiletine or 10 \lM of TTX, 
respectively reduced the INap to 0.09 ± 0.03% (n = 5, green column, 
**p= 0.002) and 0.06 ± 0.06% (n = 5, gray column, ***p= 0.001). After 
treatment, no differences were observed with WT channels. "#" indicate 
statistical differences when compared to WT channels, "*" indicate 
statistical differences when compared to Q1476R channels. 



to those in the original simulations at SS. Thus, applying solely 
Nai changes reproduced most of the APD90 decrease observed at 
SS. 

In the TNNP 2004 model, the only two ion transfer mecha- 
nisms that depend on Na; are the Na + -K + pump (INaK) and 
the Na + -Ca 2+ exchanger (INaCa) (Ten Tusscher et al., 2004). 
Their activity increase when Nai is increased. In Figure 8A, the 
two upper panels are a repeat of the first and third panels from 
the top in the middle column of Figure 7. The lower panel 
of Figure 8A is new and reports the time course of the INaK 
and INaCa currents during an AP. There was a large immediate 
increase in the duration of the first AP recorded after switching 
the model to the WR condition (Figure 8A upper panel, com- 
pare traces a and b). This was a direct effect of the increase in 
the conductance of the INap (middle panel, blue vs. black trace). 
The magnitudes of the INaK and the INaCa exchanger currents 
followed the changes in the membrane potential (compare blue 
and black traces in the lower panel). We further stimulated the 
model at 1 Hz until SS. The APD90 shortened (Figure 8A upper 
panel, trace c, orange line). While the changes in INap magnitude 
just followed the changes in membrane voltage values (orange vs. 
blue trace, middle panel), the magnitudes of both the INaK and 
the INaCa currents increased (orange vs. blue trace, lower panel). 



The increase in the INaK current was small compared to the 30% 
increase in the magnitude of the INaCa exchanger current. 

As seen in Figure 7, there was a concomitant raise in the dias- 
tolic Ca; and in the Ca; transients in tandem with the increase in 
Nai when the stimulations were applied until SS. To isolate the 
effect of the increase in Naj, we started with the WR model sta- 
bilized at 1 Hz and instantaneously changed Na; to that of the 
SS WT model (Figure 8B, black to blue traces). The model was 
then instantaneously returned to the WT model (Figure 8B, blue 
to orange traces). The changes in APD90, INaK, and INaCa were 
comparable to those in Figure 8A. This indicates that, when Nai 
was increased, the reverse mode of INaCa provided the main 
increase in the outward current that caused the shortening of the 
APD 90 . 

A possible contribution of the increase in window current in 
the effects of the mutation was tested in the model at steady- 
state at 1Hz. The diastolic membrane voltage in the WT situa- 
tion was — 85.9mV and the threshold for stimulation was 21.4 
pA/pF, whereas in the WR situation, they were — 85.5mV and 
20.8 pA/pF, respectively. Thus, there was a 0.4 mV depolariza- 
tion and a 2.8% decrease in stimulation threshold. Further, when 
the window current of the WR model was made identical to that 
in the WT model, by suppressing any changes in steady-state 
activation and inactivation curves, the diastolic membrane volt- 
age was — 85.7mV and stimulation threshold 21.3 pA/pF. Thus, 
although slightly, the window current appeared to increase cell 
excitability. 

DISCUSSION 

We biophysically characterized a novel LQT3 mutation (Q1476R) 
found in a previously symptomatic 34-year-old female with a 
classic LQT3 ECG phenotype. The Q1476 residue is a highly 
conserved amino acid in all Na + channels (Figure 2) through- 
out evolution, from humans to zebra fish. Interestingly, a second 
highly conserved Q residue (Q1475) (Figure 2) was found imme- 
diately upstream from Q1476. These two residues were located in 
the inactivation gate upstream from the IFM motif. Given their 
location, they are likely to play a role in the inactivation process. 
Moreover, the Q1489H missense substitution on SCN1A encod- 
ing the Na v l.l neuronal voltage-gated Na + channel (the Q1475 
homolog in Na v 1.5) is highly expressed in the central nervous 
system, including the retina, and is found in patients with famil- 
ial hemiplegic migraine (FHM, Figure 2B) (Vahedi et al., 2009). 
Unfortunately, no biophysical characterization has been reported 
for this mutation. We can speculate that a gain-of-function due to 
this mutation maybe the cause of FHM. 

The present paper provides support for the concept that the 
increase in the window current and the increase in the TTX sensi- 
tive persistent current due to the LQT3 mutation (Q1476R) cause 
LQTS. This is consistent with all studies indicating that muta- 
tions in SCN5A causing LQTS are related to a gain-of-function 
(Schwartz et al., 2012). 

The treatment with 50 |xM of mexiletine importantly reduced 
the INap observed on Q1476R mutant channels. This might 
explain the beneficial effects of mexiletine treatment on the 
patient. Unfortunately, the nausea due to mexiletine have 
imposed the discontinuation of the treatment. 
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FIGURE 7 | In silico simulation of heterozygous cardiomyocytes carrying 
the Q1476R mutation. (A) The three columns correspond to different 
stimulation frequencies indicated at the top of the figure. From top to bottom: 
the action potential (AP), the fast sodium current (INa), the persistent sodium 
current (INap), and cytosolic [Ca 2+ ] (Cai). Under column 0.5 Hz, traces at SS 
are superimposed: trace a: WT at 1 Hz (black), trace b: WT at 0.5 Hz 
(magenta), trace c: WR (1/2 WT + 1/2 Q1476R) at 1 Hz (orange), and trace d: 
WR at 0.5 Hz (green). Under column 1 Hz, trace a is WT at SS at 1 Hz, trace b 
is the subsequent AP after changing the model to WR, and trace c is the SS 
in WR at 1 Hz. Under column 1.5 Hz, trace a is WT at SS at 1.5 Hz, trace b is 
the subsequent AP after changing the model to WR, and trace c is the SS in 
WR at 1.5 Hz. (B) Model tests of the role of Naj in determining APD 90 . 
Instead of running the model until SS, Naj is changed instantaneously to its 
SS value. Left panel, trace a shows the SS AP in WT at 1 Hz. Trace b shows 



the subsequent AP recorded after changing Naj to the value previously 
observed at SS in the WT model at 0.5 Hz (9.84 mmol/L). Trace c shows the 
subsequent AP recorded from the model after simultaneously setting the 
model to WR and changing Naj to the value for WR at 1 Hz at SS 
(12.92 mmol/L). Trace d shows the next AP recorded after changing Naj to the 
value previously observed at SS in WR at 0.5 Hz (10.95 mmol/L). Middle 
panel, trace a shows the SS AP in WR at 1 Hz. Trace b shows the subsequent 
AP recorded after changing Naj to its value previously observed at SS in the 
WT model at 1 Hz (11.94 mmol/L). Trace c is the subsequent AP recorded 
from the model after setting the model to WT. Right panel: Trace a shows the 
SS AP in WR at 1.5 Hz. Trace b shows the subsequent AP recorded after 
changing Naj to the value previously observed at SS in the WT model at 
1.5 Hz (13.59 mmol/L). Trace c shows the subsequent AP recorded from the 
model after setting the model to WT. 



Table 2 | APD 90 , Na,, and diastolic Ca, at steady state recorded in the 
model during 1- and 1.5-Hz stimulation. 



Stimulation 


Parameter 


WT 


WR first AP 


WR 


frequency 










1 Hz 


APDgo (ms) 


454 


587 


513 




Naj (mmol/L) 


11.94 


11.94 


12.92 




Diastolic Caj (u,mol/L) 


0.109 


0.109 


0.135 




APDgo in Naj tests (ms) 


451 


582 




1.5Hz 


APDgo (ms) 


392 


471 


436 




Naj (mmol/L) 


13.59 


13.59 


14.57 




Diastolic Caj (u,mol/L) 


0.179 


0.177 


0.221 




APDgo in Naj tests (ms) 


396 


484 





As expected from its localization, the mutation impacts the 
inactivation process. The Q1476R mutation resulted in a depolar- 
izing shift of inactivation, which tended to increase the overlap of 
activation and inactivation gating (Figure 5). At voltages within 
this overlap region, Na + channels are partially but not fully 
inactivated, thus increasing the potential for persistent window 



Table 3 | Effect of bradycardia on APD 90 , Na,, and diastolic Ca, at 
steady state. 





WTat 


WTat 


WRat 


WR at 




1 Hz 


0.5 Hz 


1 Hz 


0.5 Hz 


APDgo (ms) 


454 


558 


513 


705 


Naj (mmol/L) 


11.94 


9.84 


12.92 


10.95 


Diastolic Caj (u,mol/L) 


0.109 


0.057 


0.135 


0.071 


APDgo in Naj test (ms) 




542 


515 


683 



currents (Attwell et al, 1979). At — 50 mV, the peak window 
current probability predicts that a small percentage (0.1%) of 
the mutant channels will be persistently activated. The resulting 
inward Na + current at resting membrane potentials could depo- 
larize cardiac myocytes, leading to an increase in their excitability. 
The increase in persistent current is not surprising as the muta- 
tion is localized in the inactivation gate (DIII-DIV linker). The 
effect of this mutation on the inactivation process would be a 
decrease of the probability for channels to enter an inactivated 
state. Consequently, a small amount of Na v 1.5 Q1476R channels 
will not inactivate after the rapid activation. Similar mechanisms 



www.frontiersin.org 



October 2013 | Volume 4 | Article 126 | 9 



Moreau et al. 



Sodium overload attenuate LQT3 arrhythmias 




I 1 I I I I I I 



J L 



800 



600 



O 

g 500- 

< 

400 



300 



n WT 

v WT + Q1476R 




10 



13 



15 



[Na + ], (mM) 



FIGURE 8 | Time course of INaK and INaCa and summary of SS APD 90 
and Na; values. (A,B) Analysis of the contribution of INap, INaK, and INaCa 
to APD90 changes. (A) The black traces (a) show the SS AP and currents in 
WT at 1 Hz. The blue traces (b) show the subsequent AP and currents 
recorded after introducing the heterozygous Q1476R mutation. The orange 
traces (c) were obtained after reaching the SS. (B) The black traces (a) show 
AP and currents for the WR model stabilized at 1 Hz. The blue traces (b) 
were obtained when Naj was instantanously set to Na; of the SS WT 
model. The orange traces (c) were obtained when switching to the WT 
model. (C) Summary of APD90 and Naj values at SS at three stimulation 
frequencies (0.5, 1 , and 1.5 Hz) and in the WT and WR models. The green 
arrows indicate the transitions due to bradycardia, the red arrow those due 
to tachycardia. The magenta arrows indicate the effect of passing from the 
WT to WR condition at each frequency. 



are believed to underlie the increased excitability of sensory neu- 
rons harboring inherited human pain disorder mutations that 
produce shifts in Na v 1.7 activation and inactivation of similar 
polarity and magnitude as those observed in the present study 
(Harty et al, 2006; Rush et al., 2006; Dib-Hajj et al, 2010). 

The in silico model confirmed the arrhythmogenic potential 
of the Q1476R mutation. When the model was changed to the 
heterozygous mutated paradigm, there was an immediate effect 
that prolonged APD90. Continued stimulation caused a partial 
reversion of this effect when SS was reached. At the same time, 
both Na; and Ca; increased. Figure 8C summarizes the effects 
of the mutation combined with the heart rate on APD90 and 
Na;. At SS, the accumulation of Na + ions, related both to an 
elevated frequency and to the presence of the mutation, tended 
to shorten APD90. This would in turn tend to moderate the 
LQT effect of the mutation at elevated heart rates such as under 
effort or moderate stress. The gain-of- function due to a persistent 



o 
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FIGURE 9 I Simulation of an interval variations paradigm known to 
trigger EADs and Torsades de Pointes. From SS at 0.5 Hz, a series of 40 
stimulations were applied at 666 ms intervals. Then, a combination of a 
1666 ms interval plus a couple of stimulations at 666 ms interval was 
applied twice (panel A). The action potentials produced by model for the 
last 12 s of the protocol are shown in (B). While the WT model (black line) 
shows a separate AP in response to each stimulation, the WR model 
APD90 (magenta line) was longer than the 666 ms interval, causing a 
bigeminal aspect: every second AP is smaller and shorter. After the pauses, 
EADs were present in the WR model only (arrows). 



Na + current may thus protect against arrhythmogenesis related 
to a limited increase in heart rate. However, as was seen in the 
model, an overload of Ca; stores may appear at higher rates 
and may trigger arrhythmias (Christe et al., 2008). Our data 
are consistent with published experimental and model studies 
(Pieske and Houser, 2003; Noble and Noble, 2006; Undrovinas 
et al., 2006; Soliman et al., 2012). When combined with elevated 
Ca 2+ entry under adrenergic stress, arrhythmias may appear due 
to early or late after-depolarizations. Under elevated Ca; con- 
ditions, a decrease in the L-type Ca 2+ current ought to take 
place due to the decrease in the transmembrane [Ca 2+ ] gradient 
and to the accelerated decay by calcium-dependent inactiva- 
tion. Nevertheless, its contribution to the overall effect should be 
minor. The slowly activating potassium current IKs was shown 
to be increased by elevation of Ca;, as due to a relief of KCNQ1 
inactivation by calmodulin in a Ca 2+ dependent way (Ghosh 
et al., 2006). This feature was absent in the model. Such an 
effect would enhance the frequency-dependent effect observed 
in the model, resulting into a further shortening of the AP at 
high frequencies that would be larger in the WR than in the WT 
situation. 

A combination of interval variations known to trigger EADs 
and torsades de pointes (El-Sherif et al., 1999; Viswanathan and 
Rudy, 1999; Tan et al, 2006) has been tested in the model. EADs 
appeared in the WR condition while the WT remained unaffected 
(Figure 9). This is consistent with emotional stress and exercise as 
known triggers of arrhythmias in LQT3 (Schwartz et al., 2001). 

Conversely, bradycardia causes a lower Na; and a lengthen- 
ing of APD, which are well-known features of bradycardia (Faber 
and Rudy, 2000). Changes in diastolic Ca; and Ca; transients 
are also observed. Each of these may change the magnitude of 
ionic currents and exchangers. Elevated Na; is known to cause 
an activation of the Na/K ATPase pump and an increase in the 
reverse mode Na-Ca exchanger current. Both of these processes 
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may shorten APD (Armoundas et al., 2003), as well as the increase 
in Ca 2+ -dependent inactivation of the Ca 2+ current due to ele- 
vated Ca; (Carmeliet, 2006). Our simulations showed that most 
of the APD changes were due to changes in Nai, mainly through 
a lower INaCa. The LQT effect of the mutation is likely to be 
more serious at low heart rates, as illustrated in our simulations 
where APD90 reached very critical values around 700 ms in the 
heterozygous model. Thus, the carrier of the mutation may be 
preferentially sensitive to factors triggering arrhythmias when at 
rest, which is a characteristic of LQT3 syndrome. During activ- 
ity, APD shortening, due to increased entry of Na + ions, may 
protect the carrier against arrhythmias. These results should be 
qualitatively reproduced in other models of human ventricular 
myocytes, since the basic features of the Na; homeostasis ought to 
obey to the same constraints. Furthermore, the functional conse- 
quences, as analyzed here for mutation Q1476R, may in principle 
be extended for other LQT3 mutations involving an increase in 
persistent and window currents. The intracellular Na + overload, 
as related to the gain-of-function in Na + channels, is proper 
to LQT3 and is unlikely in LQT syndromes related to potas- 
sium channel dysfunctions. This might be taken into account for 
personalized prevention and therapy of LQT3 patients. 

LIMITATIONS OF THE IN SILICO MODEL 

The basic TNNP model assumed no fraction of INap. However, 
the magnitudes of the INaCa, INaK, and Na + background cur- 
rents had been adjusted to account for the high Na; in ventricular 
human myocytes. The present modified model may thus have an 



overall excess of Na + entry. A late Na + current was added to the 
model in the WT condition, and its magnitude was limited to 
one-third of the fraction of the peak Na + current found in tsA201 
cells. This caused the APD90 in the WT configuration to be 454 ms 
vs. 325 ms in the basic TNNP model, whereas values of ~370 or 
400 ms have been reported by Li etal. (1998, 1999) in human ven- 
tricular myocytes. As a result, the effects of the mutation may have 
been overemphasized in the present simulations. 

Several ion channels have been shown to be modulated by 
intracellular Ca 2+ , either directly or through intracellular effec- 
tors (e.g., calmodulin), such modulations were not included 
in the model. Several additional channels (e.g., stretch-sensitive 
channels) found in mammalian ventricular myocytes were nei- 
ther included. A number of investigations in human ventricular 
myocytes are still needed to document their formulation. 

Lastly, this model assumes that no changes other than those 
due to the Na + channel mutation took place. There are reasons 
to doubt this, since a permanent Na + and Ca 2+ overload, as seen 
here in the WR model, are known triggers of remodeling of ion 
channel densities and excitation-contraction coupling efficiency 
in cardiac cells, which may have taken place over the long-term in 
carriers of the mutation. 
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APPENDIX 

Equations used in the TNNP model: 

WT is a variable taking three possible values: WT = 1 for a pure 
WT cell, WT = 0.5 for the heterozygous condition and WT = 0 
when all channels are mutated. 

gNaWT = WT • 14.838 



SNaQ1476R = (1 - WT) • 14.838 • (1 - 0.02) 



^infWT 



V + 86.5 

1 + e 4.4 



fynfQ1476R v + 80 . 6 

1 + e 4-3 



°- 5 '(r 



i 



138M0- 6 -e-° 1017V '+6.537-e 0 " 8016 ' 



y + 0.0005^ 



4.83 



! Na fastW T = gNaWT ■ "1WT 3 ■ ^WT ' (V ~ ^Na) 



gNapsWT 



WT - 14.838-0.0019 
fin 



1 + e -74 



fin = 3 



*Na fastQ1475R = gNaQ1476R • WQ1476R 3 • h QU76R ■ (V - £ Na ) gNapsQ1476R 



(1 - WT) • 14.838 • (1 - 0.02) • 0.0033 
fin 



1infQl476R : 



V+51.3 

1 + e -s-2 
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'(lOl.fre™ 



5 +0.022, 68-e" 
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4.83 



'NapsWT = gNapsWT • ^infWT • (V — i?Na) 
*NapsQ1476R = gNapsQ1476R ' «infQ1476R 3 ' (V — ilNa) 
*Na = *NafastWT + *NafastQ1476R + *NapersWT + NapersQli76R 
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